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Abstract 



If a point mass moves through a uniform gas cloud, at what rate does it accrete ma- 
terial? This is the question studied by Bondi, Hoyle and Lyttleton. This paper draws 
together the work performed in this area since the problem was first studied. Time 
\^q | has shown that, despite the simplifications made, Bondi, Hoyle and Lyttleton made 

quite accurate predictions for the accretion rate. Bondi-Hoyle-Lyttleton accretion 
has found application in many fields of astronomy, and these are also discussed. 
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1 Introduction 



In its purest form, Bondi-Hoyle-Lyttleton accretion concerns the supersonic 
motion of a point mass through a gas cloud. The cloud is assumed to be free 
of self-gravity, and to be uniform at infinity. Gravity focuses material behind 
the point mass, which can then accrete some of the gas. This problem has 
found applications in many areas of astronomy, and this paper is an attempt 
to address the lack of a general review of the subject. 

I start with a short summary of the original work of Bondi, Hoyle and Lyt- 
tleton, followed by a discussion of the numerical simulations performed. Some 
issues in Bondi-Hoyle-Lyttleton accretion are discussed, before a brief sum- 
mary of the fields in which the geometry has proved useful. 
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Fig. 1. Sketch of the Bondi-Hoyle-Lyttleton accretion geometry 
2 Basics 



This section is somewhat pedagogical in nature, containing a brief summary 
of the work of Bondi, Hoyle and Lyttleton. Readers familiar with the basic 
nature of Bondi-Hoyle-Lyttleton accretion may wish to skip this section. 



2. 1 The Analysis of Hoyle & Lyttleton 



Hovle and Lvttleton (1939) considered accretion by a star moving at a steady 
speed through an infinite gas cloud. The gravity of the star focuses the flow 
into a wake which it then accretes. The geometry is sketched in figure 1. 



Hovle and Lvttletonl derived the accretion rate in the following manner: Con- 
sider a streamline with impact parameter £. If this follows a ballistic orbit (it 
will if pressure effects are negligible), then we can apply conventional orbit 
theory. We have 



r9 2 



GM 



r 2 9 = (v c 



(1) 
(2) 



in the radial and polar directions respectively. Note that the second equation 
expresses the conservation of angular momentum. Setting h = £v and making 



the usual substitution u 



we may rewrite the first equation as 



de 2 



+ u 



GM 



(3) 



The general solution is u = A cos 9 + B sin 9 + C for arbitrary constants A, 
B and C. Substitution of this general solution immediately shows that C = 
GM/h 2 . The values of A and B are fixed by the boundary conditions that 



u — » (that is, r — > oo) as — > it, and that 



r = —n—— — ► — -y^ as & — »■ 7r 
d# 



These will be satisfied by 



M =^(l + CO s0)-^°sin0 (4) 



Now consider when the flow encounters the 8 = axis. As a first approxima- 
tion, the 6 velocity will go to zero at this point. The radial velocity will be v^ 
and the radius of the streamline will be given by 

1 2GM 

~r = — (5) 

Assuming that material will be accreted if it is bound to the star we have 

1 , GM 

5*- — <° 



or 



2GM ^ 

C < Chl = — — ( 6 ) 
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which defines the critical impact parameter, known as the Hoyle-Lyttleton 
radius. Material with an impact parameter smaller than this value will be 
accreted. The mass flux is therefore 



M HL = KhL^ooPoo = 3 (7) 



y oo 



which is known as the Hoyle-Lyttleton accretion rate. 



2.2 Analytic Solution 



The Hoyle-Lyttleton analysis contains n o fluid effects, which makes it ri pe for 
analytic solution. This was performed bv lBisnovatvi-Kogan et al.l ( 19791) . who 



derived the following solution for the flow field: 




GM(l + cos9) + Cv* sine 



(9) 
(10) 
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The first three equations are fairly straightforward, and follow (albeit te- 
diously) from the orbit solution given above. The equation for the density 
is rather less pleasant, and involves solving the steady state gas continuity 
equation under conditions of axial symmetry. 

Equation 4 may be rewritten into the form 



1 + e cos(6» - O 
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where e is the eccentricity of the orbit, r is the semi-latus rectum, and 9q is 
the periastron angle. These quantities may be expressed as 



e. = tan"(g|) (13) 

(14) 

r » = W < 15) 

which may be useful as an alternative form to equation 10. 

Note that these equations do not follow material down to the accretor. Accre- 
tion is assumed to occur through an infinitely thin, infinite density column on 
the 9 = axis. This is not physically consistent with the ballistic assumption, 
since it would not be possible to radiate away the thermal energy released 
as the material loses its 9 velocity. Even with a finite size for the accretion 
column, a significant trapping of thermal energy would still be expected. For 
now we shall neglect this effect. 



2. 3 The Analysis of Bondi and Hoyle 



iBondi and Hovld ( 19441 ) extended the analysis to include the accretion column 



(the wake following the point mass on the 9 = axis). We will now follow their 




Fig. 2. Sketch of the geometry for the Bondi-Hoyle analysis 

reasoning, and show that this suggests that the accretion rate could be as little 
as half the value suggested in equation 7. Figure 2 sketches the quantities we 
shall use. 

From the orbit equations, we know that material encounters the 6 = axis at 



-2 2 

2GM 



This means that the mass flux arriving in the distance r to r + dr is given by 

2nGMp 



2< dC • p 00 v c 



-dr = Adr 



(16) 



which defines A. Note that it is independent of r. The transverse momentum 
flux in the same interval is given by 



A-v 9 {0 = O) 



1 



2ns 



which is the mass flux, multiplied by the transverse velocity, divided over the 
approximate area of the wake. Applying the orbit equations once more, and 
noting that a momentum flux is the same as a pressure, we find 



Ps 



A 2GM 



2ns V r 



(17) 



as an estimate of the pressure in the wake. The longitudinal pressure force is 
therefore 



d(ns 2 P s ) = Ai 



IGM 



2 VvV 

Material will take a time of about r/v^ to fall onto the accretor from the 
point it encounters the axis. This means that we can use the accretion rate to 



estimate the mass per unit length of the wake, m, as 

GM 



msiA- 



v 3 

oo 



This makes the gravitational force per unit length 

GMmdr G 2 M 2 dr 
p rav — ~ \ — 

For accreting material, we must have r ~ GMv^. If we also assume that the 
wake is thin (s <C r) and roughly conical (ds/s ~ dr/r), then taking the ratio 
of the pressure and gravitational forces, we find that pressure force is much 
less than the gravitational force. We can therefore neglect the gas pressure in 
the wake. 

The mass per unit length of the wake, m, was introduced above. If we assume 
the mean velocity in the wake is v, we can write two conservation laws, for 
mass and momentum: 



— {mv) = A (19) 

dr 

d , 9 . GMm , . 

— (mv 2 )=Av oc — (20) 

dr r z 

Recall that Ar^ is the momentum supply into the wake, since r = v^ on axis 
for all streamlines. We can declutter these equations by introducing dimen- 
sionless variables for m, r and v. 



(21) 

GM , . 

r = —X (22) 

"to 

AGM . . 

m— — ^ — /i (23) 
v 

Note that x = 2 corresponds to material arriving from the streamline char- 
acterised by £ HL . Substituting these definitions into equations 19 and 20, we 
obtain 



A M = i (24) 

A(^2 ) = 1 _^-2 (25) 

We shall now analyse the behaviour of these equations. 
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We can integrate equation 24 to yield 

ixv = x ~ a (26) 

for some constant a. Since \x is a scaled mass (and hence always positive), we 
see that the scaled velocity [y) changes sign when x — ol. That is, a is the 
stagnation point. Material for x < a wm accrete, so knowing a will tell us 
the accretion rate (since the accretion rate will be Aro where r$ is the value 
of r corresponding to a). By writing jj,u 2 — /iv • u, we can use equation 26 to 
rewrite equation 25 as 

dx X-ot X 

This has not obviously improved matters, but we can now study the general 
behaviour of the function, without trying to solve it. First we need some 
boundary conditions. These are as follows: 



v — *• 1 as x —* °° 




Which means that v — 


* v^ at large radii 


v = at x — a 




The stagnation point 





— > Everywhere 
dx 

The velocity must be a monotonic function. This is physically reasonable, 
if we are to avoid unusual flow patterns 

The first two conditions can be satisfied for any value of a. Fortunately, the 
third implies as restriction. The next set of manipulations may seem a little 
obscure at first, but they do lead in the desired direction. 

Substitute £ = a~ 1 x- Equation 27 then reads 

d£ £-1 a? [ ) 

Now, suppose the derivative is zero. This leads to the condition 

v 2 -v + -L(x-l) = 
or, one application of the quadratic roots formula later: 



Wi-^*- 1 ' (29) 
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Fig. 3. Curves where 4j = for a < 1. In the regions marked 'a,' the derivative is 
greater than zero. It is less than zero in the 'b' region 

Since v ultimately represents a physical quantity (the velocity), it's obviously 
desirable that it remain real. We therefore need to look at when the discrimi- 
nant can become zero. This is another quadratic equation, leading to 



i = - (i ± vr^) 



which means that something must happen when a 
'something,' it is best to plot equation 29. 



1. To determine this 



Figures 3 and 4 demarcate the regions where ^| changes sign, as dictated 
by equation 29. These are not possible solutions for v. However, any suitable 
solution for v must remain within the region marked 'a' if it is to remain 
monotonic and increasing. This is only possible when a > 1. Unwinding our 
rescaled variables, we see that an a value of unity puts the stagnation point 
halfway between the accretor and the original value of iHovle and Lvttleton . 
This in turn implies a minimum accretion rate of 0.5M HL . 



Again, I would like to remind the reader that the flow has been assumed to 
remain isothermal with negligible gas pressure throughout this discussion. This 
assumption is likely to be violated in the wake, where densities will be high 
and radiative heat loss inefficient. At the very least, thermal effects should be 
important close to the stagnation point in the wake. Ifioredd ( 20001 1 details an 
analysis similar to that gi ven above, but with a pressure term included. The 
value of a (which iHoredtl calls xo) was found to lie between 0.6 and 3.5 for 
flows which were supersonic at infinity and subject to Newtonian physics (the 
polytropic and adiabatic indices were also free parameters in this analysis). 




Fig. 4. Curves where 4? = for a > 1. In the region marked 'a,' the derivative is 
greater than zero. It is less than zero in the 'b' regions 



Will the flow be stable? Bondi an d Hoyle) asse rted that if a > 2 (note that 
a = 2 gives the solution of lHovle and Lvttleton ). then the wake would become 
unstabl e to pe r turba tions which preserve axial symmetry. However, later anal- 
ysis bv ICowid ( 1977 ) suggested that the wake should be unstable, regardless 
of the value of a. Subsequent numerical simulations and analytic work have 
shown that Bondi-Hoyle-Lyttleton flow is far from stable, and we will discuss 
the subject in section 4.2. 



2.4 Connection to Bondi Accretion 



Bondil ( 19521 ) studied s pherically symmet ric accretion onto a point mass. The 
analysis shows (see e.g. lFrank et al.l ( 20021 )) that a Bondi radius may be denned 



as 



rii 



GM 
c 2 s (r B ) 



(30) 



Flow outside this radius is subsonic, and the density is almost uniform. Within 
it, the gas becomes supersonic and move s towar ds a freefall solution. The 
similarities between equations 6 and 30 led iBondil to propose an interpolation 
formula: 

2ttG 2 M 2 Poo 



M 



(cl + ^) 3 / 2 



(31) 



This is often known as the Bondi-H o yle ac cretion rate. On the basis of their 
numerical calculations, IShima et al.l (|l985[ ) suggest that equation 31 should 



acquire an extra factor of two, to become 

M - ^ G2M2 PQQ (32) 

which then matches the original Hoyle-Lyttleton rate as the sound speed be- 
comes insignificant. The corresponding £ BH is formed by analogy with equa- 
tion 7. 

Nomenclature in this field can be a little confused. When papers refer to 
'Bondi-Hoyle accretion rates,' they may mean equation 7, 31 or 32. In this 
review, I shall refer to pressure-free flow as 'Hoyle-Lyttleton' accretion and 
use M HL and ( RL . When there is gas pressure, I will talk about 'Bondi-Hoyle 
accretion' and use M BH and £ BH , in the sense defined by equation 32. 1 shall use 
'Bondi-Hoyle-Lyttleton' accretion to refer to the problem in general terms. 



3 Numerical Simulations 



In the previous section, I outlined the basic theory behind Bondi-Hoyle- 
Lyttleton accretion. This lead to elegant predictions for the accretion rate, 
as given by equations 7 and 32. However, reaching these equations required a 
lot of simplifying assumptions, so necessitating further investigation. The in- 
tractability of the equations of fluid dynamics requires a numerical approach 
to the problem. 

In a break with tradition, I shall start this section with the answer, and then 
give more detailed citations to examples. 



3. 1 Summary 



Do the equations of iBisnovatvi-Kogan et al.l provide a good description of 



Bondi-Hoyle-Lyttleton flow? The answer is 'No.' 

In the absence of other effects, three numbers parameterise Bondi-Hoyle- 
Lyttleton flow: 

• The Mach number, M. 

• The size of the accretor, in units of £ HL 

• The 7 value of the gas 

Figure 5 shows sample density contours for a flow with Ai = 1.4, an accretor 
radius of 0.1£ HL an d 7 = 5/3. This particular simulation was axisymmetric. A 
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Fig. 5. Density contours for a sample Bondi-Hoyle-Lyttleton simulation. The flow 
had M = 1.4, an accretor radius of 0.1£ HL and the equation of state was adiabatic 
with 7 = 5/3. The contours are logarithmically spaced over a decade of density. The 
dotted line indicates £ HL . The flow is incident from the left 
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Fig. 6. Velocity field corresponding to the densities shown in figure 5. The approxi- 
mate position of the bow shock is marked with a dotted line 

bow shock has formed on the upstream side. The corresponding velocity field 
is plotted in figure 6. Downstream of the shock, material flows almost radially 
onto the accretor, in marked contrast to the analytic solution of equations 8 
to 11. 

But what of the accretion rate? Figure 7 shows the accretion rates obtained 
for three simulations. Although the dimensionless parameters were kept the 
same, the physical scales and grid resolution varied. Figures 5 and 6 were taken 
from run 2. The accretion rates achieved are quite close to the value of M HL 
predicted for the flow (this value is substantially larger than the corresponding 
M BH ). Despite the simplifications made, the work of Bondi, Hoyle and Lyttle- 
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Fig. 7. Accretion rates for plain Bondi-Hoyle-Lyttleton flow. The crossing time 
corresponds to C HL 

ton has been largely vindicated. In the remainder of the section, I shall cite 
places in the literature where further simulations of Bondi-Hoyle-Lyttleton 
flow may be found. 



3.2 Examples in the Literature 



Huntl computed nume rica l solut ions of Bondi-Hoyle-Lyttleton flow in two pa- 
pers written in Il97ll and Il979t The accretion rate suggested by equation 32 
agree d well with that observed, despite the flow pattern being rather different. 



Huntl studied flows which were not very supersonic and were non-isothermal. 



A bow shock formed upstream of the accretor. Upstream of the shock, the flow 
pattern was very close to the original ballistic approximation. Downstream, 
the gas flowed almost radially towards the point mass. A s ummary of early ca l- 
culations of Bondi-Hoyle-Lyttleton flow may be found in lShima et al.l (1985). 
The calculations in this paper are in broad agreement with earlier work, but 
some differences are noted and attributed to resolution differences. 



More rece ntly, a s eries of calculations i n three dimens i ons have been per- 
formed bv IRuffertl in a series of p apers ([Ruffertl . Il994at iRuffert and Arnett , 
19941 IRuffertl . Il994hl Il995l Il99ffl . This series of papers used a code based 



on nested grid s, to permit high resolution at minimal computational cost. 
IRuffertl ()1994a ) details the code, and presents simulations of Bondi accretion 
( where the accretor is s tation ary) . Bondi-Hoyle-Lyttleton flow was considered 
in JRuffert and Arnettf (|l994J ). The flow of gas with M = 3 and 7 = 5/3 past 
an accretor of varying sizes (0.01 < r/( BU < 10) was studied. For accretors 
substantially smaller than £ BH , the accretion rates obtained were in broad 
agreement with theoretical predictions. The flow was found to have quiescent 
and active phases, with smaller accretors giving larger fluctuations. However, 
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these fluctuations were far less violent t han th e 'flip-flop' instability observed 
in 2D simulations (see below). iRuffertl (1994b) extended these simulations to 



cover a range of Mach numbers, finding that higher Mach numbers tended 
to give lower accretio n rates (down to the original interpolation form ula of 
equation 31). IRuffertl studied the flow of a gas with 7 = 4/3 in the Il995l 



paper, finding accretion rates comparable with the theoretical results. Small 
accretors and fast flows were required before a ny instab i lities appeared in the 
flow. Nearly isothermal flow was considered in IRuffertl (J1996). The accretion 



rates were slightly higher than the theoretical values (except for the smaller 
accretors), and the shock moved back to become a tail shock. The oscillations 
in the flow were less violent still. 

The reason for the formation of the bow shock is straightforward - the ris- 
ing pressure in the flow. As shown by equation 11, the flow is compressed as 
it approaches the accretor. This compression will increase the internal pres- 
sure of the flow, eventually causing a significant disruption. At this point, the 
shock will form. This interpretation is consistent with the behaviour observed 
in simulations, where decreasing 7 moves the shock back towards the accre- 
tor. However, the precise location of the shock does no t seem to be a strong 



function of the Mach number (cf the papers of IRuffertl ) 



4 Issues in Bondi-Hoyle-Lyttleton Flow 



In this section I shall discuss some issues relating to Bondi-Hoyle-Lyttleton 
flow which are of particular interest. 



4.1 The Drag Force 



The simple idealisation of Bondi-Hoyle-Lyttleton flow cannot persist for long. 
The accretor is not only increasing its mass - it is accumulating momentum 
as well. Eventually, it should be accelerated to being co-moving with the gas 
flow. A full calculation is not straightforward, but dimensional considerations 
suggest 

M Voo ~ M Voo = F drag (33) 

Please note, that the ~ in this equation is very approximate. However, equa- 
tion 33 suggests that the accreting body will be brought to rest with respect 
to the flow on the mass doubling timescale. This is obviously a problem if the 
accretor is to change its mass appreciably. As we shall see later, this has lead to 
most research into Bondi-Hoyle-Lyttleton accretion being concentrated into 
the study of binaries. In such cases, the momentum difference can be 'paid' 
by a change in orbit. 
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The drag does not originate as a form of 'wind' resistance pressing directly on 
the accretor. This is for two reasons 



• The momentum deposited by the accretion column will be far larger 

• Mathematically, the accretor is a point anyway 

Instead, the drag arises from the gravitational focusing of material behind the 
accretor. Since more material is present on the downstream side, the gravita- 
tional attraction of th e downstream side is larger, l and exerts a drag force. 
Chandrasekharl ( 19431 ) was the first to consider this problem - called 'dynam- 
ical friction' - for a collisionless fl uid. An extended (and more recent) dis- 
cussio n of the problem is given by iBinnev and Tremaind ( 19871 ) . iDokuchaev 
(J1964J) discussed t he problem for a gaseous me dium. The matter of drag is 
also mentioned bv lRuderman and Spieeell ( 19711 ). who propose 



drag 



M Voo \n 



(34) 



where b min and 6 max are cut-off radii for the gravitational force. lYabushita 



119781 ) suggested that a suitable value for the outer cut-off for a flow with 
pressure would be the point where the pressure in the wake became equal to 
the background pressure (the inner cut-off radius is usually taken to be the 
radius of the accretor itself). 



Values for the drag force given bv IShima et al.l ( 1985T ); IShankar et al.l (1993); 
Rufiert and Arnettl (J19941 ) suggest that the drag force is no more than a factor 
of ten larger than the crude estimate of equation 33. The precise drag value has 
a tendency to fluctuate anyway - Bondi-Hoyle-Lyttleton flow is not stable. 



4-2 Flow Stability 



Even in the axisymmetric case, there is no particular reason to believe that 
Bondi-Hoyle-Lyttleton flow should be stable. The binding energy test of equa- 
tion 6 is made for gas flowing away from the accretor. If this material is going 
to be accreted, it needs to turn around somehow and fall towards the point 
mass. This m ust happen in some sort of accretion column, of the ty pe first 
consid ered bv iBondi and Hovld ( 19441 ) . As noted above, the work of ICowid 
( 1977T ) found that this wake should be unstable. A 'shock cone' must surround 
the wak e, in which th e flow loses its 9 velocity before it encounters the axis (see 
also [Wolfsonl (|l977al )b[)). The high densities expected for the wake mean that 
this shock is likely to heat the gas. Gas pressure could then be expected to 



1 This is ignoring the mathematical impossibility of an infinite, uniform medium 
for the unperturbed case 



14 



drive oscillations close to the stagnation poi nt. Bondi-Hoyle-L vttleton flow 
around small accretors has been studied by JKoide et al.l (J19911) . This paper 
notes that that the 'accreting body is so small that a part of the accreting gas 
sometimes misses the target object and flows towards the upstream as a jet.' 
This is obviously a rather unstable condition, and leads to the accretion flow 
'sloshing' back and forth around the point mass. The accretion rate fluctuates 
too, although the time averaged rate is still close to the Bondi-Hoyle value. 



When th e condition of axisymmetry is relaxed, even more i nstabilities become 



possib le. iMatsuda et al.l (|1987[ ): iFrvxell and Taaml (|1988| ): iTaam and Frvxelll 
( 19881 1 performed two dimensional simulations of the Bondi-Hoyle-Lyttleton 
geometry, with the condition of axisymmetry relaxed, and a density and/or 
velocity gradient imposed on the upstream flow. All three papers found that 
a 'flip-flop' instability resulted, with the wake osci llating back and forth in a 
manner reminiscent of a von Karman vortex street. IMatsuda et al.l ([19911 ) sug- 
gested that the instability was intrinsic to the accretion fl ow, since it was found 
to develop even under the conditions first considered by Hoyle and Lyttleton 
. Another detailed study of the ' flip-flo p' instability for the 2D case for 



(1939). 



isothermal gas is that of IShima et al.l ( 1998J ). The code used was specifically 
designed to conserve angular momentum and to permit very high resolution 
in the inner portions of the grid. They suggest that some of the resolution de- 
pendence of the instability found by earlier work was due to the use of codes 
which conserved linear momentum, and caution against the results obtained 
by such codes. 

Such a spe ctacular insta bility naturally prompted a n intense the oretical in- 
vestigation. |SokerJ (J1990) extended the earlier work of lCowid ( 19771 ) to include 
tangential oscillations. The analysis is based on the assumption of a pressure 
free flow, and the expressions derived for the tangential behavio ur also re- 
quire the flow to be 2D and planar. The radial instability noted bv ICowid was 
found to be independent of the incoming material. Any radial oscillation in 
the wake would grow, although the growth timescale was much longer than 
the oscillation timescales. The tangential mod es (cor responding to the 'flip- 
flop' instability) behaved in a similar manne r. ISokerl also predicted that the 
instability should be milder in the 3D case. In lSokerl ([1991I ). a numerical study 
of the coupling between the radial and tangential oscillations was made. The 
radial modes, corresponding to large density and velocity fluctuations were ex- 
cited far beyond the linear regime, while the tangential oscillations remained 
linear. While the mass accretion rate showed corresponding fluctuations, the 
time averaged accr etion rate was simi lar (although smaller) than the predic- 
tion of equation 7. ILivio et al.l (|1991[ ) added a simple analysis of the shock 
cone surrounding the wake seen in numerical simulations. Instabilities were 
found in both the planar 2D and full 3D cases, although the authors not e that 
the in stability should be milder in the 3D case (a point also made by ISoker 
(GHU). 
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The major weakness of all simulations of the 'flip-flop' instability mentioned so 
far is that they fundamentally change the geometry of the problem. In order to 
simulate non-axisymmetric flow in two dimensions, the flow has to be assumed 
to be planar. This changes the shape of the accretor from a sphere to a cylinder. 
The equations of fluid dynamics are non-linear, and are notorious for thei r 
resolute refusal to yield to a proof of solution uniqueness ( Feffermanl |2000). 
There is therefore no particular reason to expect the 2D plana r simula tions to 
be characteristic of the true solution in 3D. The simulations of iRuffertl suggest 
that the 'flip-flop' instability is an artifact of 2D planar flow. 



Foglizzo and Ruffertl ( 19971 Il999j ) att empted to model the instabilities ob- 
served in the earlier numerical work of lRuffertl . The first of these papers con- 
structs stationary models, while the second contains a stability analysis. The 
origin of the instability was the bow shock generally seen in numerical simu- 
lations. This produces entropy gradients in the flo w, which allows Rayle igh- 
Taylor and Kelvin-Helmholtz instabilities to grow. JFoglizzo and Ruffertl con- 
cluded that the instability should be stronger if 

• The shock is detached from the accretor (as is the case for higher 7 values) 

• The flow has a higher Mach number 

• The accretor is smaller 



They found that the instability should be non-axisymmetri c , and start at 
around 9 = tc/2 and close to the accretor. iFoglizzo and Tagged ( 20001 ) describe 
the instability as 'entropic-acoustic,' where entropy perturbations introduced 
by a shock propagate back to the shock via sound waves. These then trigger 
new entropy perturbations. 



4-3 Non-uniform Boundary 



Suppose the conditions at infinity are not uniform, but instead a density 
and/or velocity gradient is present. This means that the flow within a cylinder 
of radius £ HL possesses angular momentum about the accretor. How much of 
this reaches the accretor? 



Early calculations (IDodd and McCrea . 1952t llllarionov and Sunvaevl . Il975t 
Shapiro and Lightmanl . Il976t IWanel Il981) suggested that most of this an- 
gular momentum would accrete. However. iDavies and Pringld (1980) pointed 
out that only material which had lost most of its angular momentum would be 
able to settle onto a small accretor. They developed a simple analytic model 
for small density and velocity gradients (if the gradients become large, then the 
flow ceases to behay e in the manner described by Bondi, Hoyle and Lyttleton). 
Davies and Pringld found that the mass accretion rate should be unaffected, 
and there should be no accretion of angular momentum. 
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Such confusion calls for numerical work. iRuffert and Anzerl ( 19951 ) presented 



a sample 3D simulation of accretion under such conditions. T hey found that 

about 70% of the angular m omentum available (as calculated bv lShapiro and Lightman 
(11976r0 would be accreted. iRuffertl e xtende d this calculation with papers in 
Il997l (covering velocity gradients) and ll999l (studying density gradients). The 



mass accretion was not affected much, while the angular momentum accretion 



rate v aried between 0% and 70% of the value suggested bv lShapiro and Lightman 



( 19761 ). Smaller accretors gave less stable flow, but none were as viole ntly un- 



stable as the 'flip-flop' instability observed in 2D planar simulations. IRuffert 
noted that a very small accretor will not be able to accrete all of the angular 
momentum within ( RL , but was unable to test such a case, due to the vast 
computational load involved in such a simulation. 



4-4 Radiation Pressure 



Where there is accretion, there will be an accretion luminosity. 2 Radiative 
feedback has the potential to alter the Bondi-Hoyle-Lyttleton flow, and a 
number of workers have studied this. 

Most work has concentrated on the problem of radiative feedback in X-ray 
binaries. This is a fairly straightforward application, since the constancy of 
the Thompson cross section makes the transfer problem intrinsically grey. 



iBlondin et al.l ( 19901 ) simulated a compact object accreting an O star wind 



(forming an X-ray binary). Radiative heating was included (it was relevant to 
the radiation-driven wind), but the radiation force of the X-rays was neglected. 
The gas was assumed to be optically thin and in ionisation equilibrium. The 
wake was found to be unstable, oscillating back and forth. Despite these, the 
accretion rates were broadly consistent with that expected from a Bondi-Hoyle 
type analysis (some modifications were necessary to allow for t he geometry 



of the binary). A later study of the same problem was made bv iTaam et al 



( 199ll ). This work included the effect of radiation pressure, but the flow was 



still assumed to be optically thin to electron scattering. Radiation pressure 
was then negligible except in the wake. This paper contained short section 
considering the accretion of an optically thin gas subject to radiation pressure. 
This led to the prediction that the flow would be unstable to oscillations if 
the accretion rate exceeded one third of the Eddington Limit. 



Kiev et al.l ( 19951 ) were interested in the application of Bondi-Hoyle-Lyttleton 



accretion following a nova explosion in a binary. The radiation field was sim- 
ulated using flux limited diffusion, and analytic approximations to opacity 



Barring certain cases of finely tuned accretion onto a black hole 
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values. Their simulations had quite a complicated model for the accretor, in- 
cluding an envelope. Radiation pressure was found to be critical to simulating 
the flow (heating was included in all calculations, but the radiation pressure 
was omitted from some). For hot, optically thick flow, including radiation 
pressure made the flow subsonic and substantially reduced the drag. However, 
the accretion rate was found to be low (much less than the predictions of 
Bondi-Hoyle theory) in all cases. 



4-5 Relativity 



If the accreting object is a neutron star or a black hole, relativistic effects are 
likely to become important. However, relativistic hydrodynamics is generally 
recognised to be a non-trivial problem, and relatively little work has been done 
on Bondi-Hoyle-Lyttleton accretion for relativistic flows. 



An early study by IPetrich et al.l (119891) found broad agre ement between the 



Newtonian and relativistic cases. iFont and Ibanea (|1998b|) performed axisym- 



metric calculations in a Schwarzschild metric. The passage of time gave more 
pow erful computer s, enabling simulations to be run with higher resolution 



than Petrich et al.. Most of the accretion rates found were similar to the New- 



tonian estimates, but so me were an order of magn itude or so higher. No signs 
of instability were found. IFont and Ibanea ( 1998al ) relaxed the assumption of 



axisymmetry, but still found the flow to be steady. However, the authors note 
that they were unable to push the simulation parameters very far. 



Font et al.l f 1999i ) simulated Bondi-Hoyle-Lyttleton flow onto a rotating (Kerr) 



black hole. Two forms of the metric were used, to differentiate between nu- 
merical and physical effects. They found that rotational effects were confined 
to a region close to the hole. Their flows remained steady. 



5 Applications 



The Bondi-Hoyle-Lyttleton scenario has been applied to a variety of problems 
in the years since its introduction. I will now discuss some examples of these. 
Each of these problems could be the subject of a review article by themselves 
(and often have been), so I am not able to discuss the issues involved in great 
depth. However, I hope that this will provide a reasonable sample of the wide 
range of areas where the Bondi-Hoyle-Lyttleton model has proven useful. 



5.1 Binary Systems 



The problem of accretion in a binary system seems to be the most popular 
application of the Bondi-Hoyle-Lyttleton analysis. As noted above, the in- 
evitable drag force simply causes the orbit to change. The gas supply can 
either be from a stellar wind, or from common envelope (CE) evolution. 



Wind accretion seems to be one of the most popular applications of the Bondi- 
Hoyle-Lyttleton geometry. However, there are a number of potential compli- 
cations. For the Bondi-Hoyle-Lyttleton solution to be valid, accreti on must be 
driven by a wind, rather than by Roche lobe overflow. The work of lPetterson 
(1978|) warns that the presence of Roche lobe overflow will substantially com- 
plicate matters, and that allowance must be made for the possibility when 
comparing theory with o bservations. The orbital mo ti on of the binary can 
also cause problems - see iTheuns and Jorissenl ( 19931 ) ; iTheuns et al.l ( 19961) 
for a discussion. These papers studied accretion rate in a binary, where the 
wind speed was comparable to the orbital velocity. Consequently, the flow 
pattern w as substantially differ ent from that envisaged by Bondi, Hoyle and 
Lyttleton. ITheuns et al.l (1996) found that the accretion rate for a binary is 
decreased by a factor of about ten compared to the prediction of equation 32. 
They attribute this difference to the disrupt i ng ef fect of the orbital motion. 
Moving to higher mass ratios, IStruck et al.l ( 2004T ) modelled wind accretion 
onto substellar companions of Mira variables. They found that the mean ac- 
cretion rates were generally similar to those predicted by Bondi, Hoyle and 
Lyttleton. However, the flow was highly variable - both due to instabilities 
in the accretion flow and the intrinsic variability of the star. The ultimate 
'wind' is that produced by an explosion. In this case, the very high velocities 
involve d tend to make Cn u comparable to the size of the accreting body. The 
work of Kiev et al.l ( 19951 ) has already been m entioned. However, the problem 
had been studied before - iMacDonaldl ( 19801 ) estimated accretion rates using 
modified a Bondi-Hoyle-Lyttleton formula. 



Jackson! ( 1975J) used the predictions of Bondi-Hoyle-Lyttleton the ory to derive 



the sy s tem parameter s of Ce n X- 3. Similar work was p erformed bv lEadie et al. 



( 1975J) , IPounds et al.l ( 1975J) and lLamers et al.l ( 19761) . By consider ing the ac- 
cretion of a high velocity wind by neutron stars. IPfahl et al.l ( 20021 ) concluded 
that most of the low luminosity, hard X-ray sources known in our Galaxy could 
be powered by such systems. Modified Bondi-Hoyle-Lyttleton accretion has 
also been used to study cases where a giant star's wind is being a c creted 
by a main sequence star. Some examples a re the work of IChapmanl (|198ll) 
and IChe-Bohnenstengel and Reimera ( 19861 ). The Bondi-Hoyle-Lyttleton ge- 
ometry is a useful first approximation to wind accretion in binary systems. 
However, unless the wind speed is much greater than the orbital velocity, the 
accretion rates can deviate significantly from the simple predictions. 
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At the extreme end of the mass ratio scale, Bondi-Hoyle-Lyttleton accretion 
has eve n been used to e s timat e accretion rates onto a planet embedded in 
a disc ( Nelson and Bend . I2003T ). Although the situation simu lated was not 



entirely appropriate to the original analysis, iNelson and Benzl point out that 
it represents a maximum possible accretion rate. This rate turns out to be 
extremely high, showing the need for higher resolution simulations of planetary 
accretion flows. 

Common Envelope evolution occurs when two stellar cores become embed- 
ded in a large gas envelope. Such an envelope is typically produced when one 
of the members of the binary system swells as it leaves the main sequence. 
For a more detailed d iscussion of Common Envelope evolution itself see, e.g. 
Iben and Livid (|1993J ). In such cases, accretion rates are critical for determin- 



ing the detailed evolution of the system. In computing accretion rates, modified 



Bondi -Hoyle-Lyttleton formul ae are often used e.g. iTaam and Bodenheimer 
( 19891 1 - see also the review bv ITaam and Sandquisd ( 20001 ). 



5.2 Protostellar Clusters 



Bondi-Hoyle-Lyttleton flow is also likely to be applicable to regions of star 
formation. Although single stars will be stopped by the drag force, real stars 
are generally thought to form in clusters. Protostars and gas are trapped inside 
a gravitational well, and orbit within it. The drag will then simply cause a 
change in orbit - as is the case for X-ray binaries. Indeed, the approximation is 
likely to be better for protostars in a protocluster. This is because the orbital 
motion of the protostars is the 'source' of the wind, rather than a wind from 
a companion. As a result, the geometry is simpler (since the orbital motion 
does not have to be added to the wind velocity). Furthermore, non-inertial 
forces (coriolis and centrifugal) are likely to be far less significant. 



Bonnell et al.l ([20011 ) performed a thorough study of accretion in a protc 



o- 

cluster. They simulated the evolution of a gas cloud containing many small 
point masses (representing protostars). The point masses gr ew by accret 



ing the gas. When the gas dominated the mass of the cluster, iBonnell et al 
found that the accretion was best described by tidal lobe overflow (examin- 
ing whether material was bound to the cluster or the star). However, as the 
mass in stars grew, Bondi-Hoyle-Lyttleton accretion became the more sig- 
nificant mechanism. The transition occurred first for the most massive stars 
which ha d sunk into the cluster c ore. However, massive stars are very lumi- 
nous and lEdear and Clarke! ( 20041 ) showed that radiative feedback can disrupt 



the Bondi-Hoyle-Lyttleton flow once stellar masses exceed 10 M Q (this is ob- 
viously dependent on the prevailing conditions in the protocluster). 
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Unfortunately, direct observations of this process are not available. Protoclus- 
ters contain large quantities of dusty gas, which greatly obscure regions of 
interest. Furthermore, the expected luminosities are lower, and the emission 
wavelengths less distinctive than those of X-ray binaries. 



5.3 Galaxy Clusters 



Galaxies orbiting in a cluster are another candidate for Bondi-Hoyle-Lyttleton 
accretion. One immediate complication is the high temperature of the inter- 
galactic medium (IGM). The IGM is typically hot eno ugh to emit X-rays, and 
hence the galactic motions will usually be subsonic. iRuderman and Spiegell 
(1971) suggested that the IGM might be heated (at least in part) by the 
accretion shocks inherent to Bondi-Hoyle-Lyttleton accretion. Galaxies are 
also rather porou s objec ts, and contain their own gas. In a study of M86, 
Rangaraian et al.l ( 19951 ) concluded that the 'plume' observed was proba- 
bly the result of ra m-pressure stripping of material from the galaxy itself. 
Stevens et al.l (J1999J) simulated galaxies under such conditions, and concluded 
that "the ram-pressure stripped tail will usually be the most visible feature." 
This paper also contains a list of observed wakes. 



De Young et al.l ( 19801 ) observed M87, and found evidence for subs onic Bondi 



Hoyle-Lyttleton flow. However, higher resolution observations by Owen et al. 



( 20001 ) suggest that this simple picture is not sufficient. In particular, the active 
nucleus of M87 drives an o utflow. A filament h as been observed trailing Abell 
1795 both in the optical ( Cowie et al. . 19831 ) and in X-rays ( Fabian et al. . 
2001). It has been proposed ( Sakellio u et all Il996l ) that this filament is an 
accretion wake, but lFabian et al.l ((2001) note that the gas cooling ti mes aren't 



quite right for this simple approximation to be completely valid. ISakelliou 
(2000) constructed a simple theoretical model of a Bondi-Hoyle-Lyttleton 
wake behind a galaxy. Wakes were expected to extend for up to 20 kpc. The 
wakes would form behind slow moving, massive galaxies in low temperature 
clusters. 



5.4 Other Applications 



The original application of Bondi-Ho yle-Lyttleton accretion w as to the flow of 
the interstellar medium past the Sun. ISikivie and Wic k (2002) apply a similar 
analysis to the flow of dark matter past the Sun. They suggest that annual 
variations in WIMP detections may be partially attributable to the focusing 
of the flow by the Sun. 



Bondi-Hoyle-Lyttleton accretion was invoked bv iKamp and Paunzenl ( 20021 ) 
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to explain the unusual chemical abundances of A-Bootis type stars. These stars 
have metal abundances typ ical of the interstellar medium (that is, metal- 
poor). iKamp and Paunzenl suggest that radiation pressure on dust grains 
(which contain most of the metals) prevents the accretion of the heavier el- 
ements, while gas accretes in a Bondi-Hoyle-Lyttleton fashion. Moving to 
higher energies, accretion onto neutron stars moving through gas clouds has 
been proposed as a m echanism for producing X-ray sources in the Galax y 
( Ostriker et all Il970f ) and in globular clusters ( Pfahl and RappaportL 12001 ). 



However with neutron stars, magnet ic fields can cause significant complica- 



tions - see, e.g. iToropina et al.l (J200l|). 



Maeda et al.l ( 20021 ) studied the central portions of our Galaxy with Chandra. 



Finding evidence for recent activity, they suggest that this could have been 
powered by the central black hole accreting material from an expanding su- 
pernova shock. This would be a transient example of Bondi-Hoyle-Lyttleton 
accretion. The potential luminosity from this is rather high (comparable with 
the Eddington Limit). However, there is a complication due to the thermal 
pressure of the ambient gas, which could reduce the inferred luminosity sub- 
stantially. 



6 Summary 



In their original analyses, Bondi, Hoyle and Lyttleton made many simplifi- 
cations. Despite these, the broad picture they present seems to be correct. 
Numerical studies have been made of the purely hydrodynamic problem, and 
of cases where extra physical processes are relevant. Bondi-Hoyle-Lyttleton 
accretion has also been used to explain phenomena in a variety of astronomical 
contexts. 

Of course, the original Bondi-Hoyle-Lyttleton results cannot be applied with- 
out some thought. Numerical studies have shown that the flow pattern is 
more complicated than that originally envisaged. Meanwhile, real systems are 
always more complicated than theoretical ones. Bondi-Hoyle-Lyttleton accre- 
tion should be regarded as a reference model - it is unlikely to explain any 
system in detail, but it can serve as a useful basis for classifying behaviours. 
It can be applied as a test model on systems of all scales - from binary stars 
up to galaxies in clusters. 

There are many future avenues for research. As well as improving simulations 
of accretion in binaries, studies need to be made on Bondi-Hoyle-Lyttleton 
accretion for flows where the accretor is orders of magnitude smaller than the 
accretion radius. Previous numerical work has hinted at lurking instabilities, 
but sufficient computing resources have yet to be brought to bear. Extra phys- 
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ical processes are also candidates for inclusion, especially radiative feedback. 
Radiative feedback is likely to be relevant in the context of both X-ray binaries 
and star formation simulations, and could alter the flow pattern significantly. 
On the observational side, better observations of systems previously modelled 
with a simple Bondi-Hoyle-Lyttleton analysis will show deviations, which can 
then be used to enhance our understanding of them. 

Time has not eroded the usefulness of the Bondi-Hoyle-Lyttleton accretion 
geometry. It provides a simple framework for examining and refining theory 
and observation. 
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